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ABSTRACT 

Rapid spectral changes in the hard X-ray on a time scale down to ~ 0.1 s are studied by applying 
"shot analysis" technique to the Suzaku observations of the black hole binary Cygnus X-1, performed 
on 2008 April 18 during the low/hard state. We successfully obtained the shot profiles covering 10- 
200 keV with the Suzaku HXD-PIN and HXD-GSO detector. It is notable that the 100-200 keV 
shot profile is acquired for the first time owing to the HXD-GSO detector. The intensity changes in 
a time-symmetric way, though the hardness does in a time-asymmetric way. When the shot-phase- 
resolved spectra are quantified with the Compton model, the Compton ^/-parameter and the electron 
temperature are found to decrease gradually through the rising phase of the shot, while the optical 
depth appears to increase. All the parameters return to their time-averaged values immediately within 
0.1 s past the shot peak. We have not only confirmed this feature previously found in energies below 
~ 60 keV, but also found that the spectral change is more prominent in energies above ~ 100 keV, 
implying the existence of some instant mechanism for direct entropy production. We discuss possible 
interpretations on the rapid spectral changes in the hard X-ray band. 
Subject headings: accretion, accretion disks — X-rays: binaries — X-rays: individual (Cyg X-1) 



1. INTRODUCTION 

Starting with the first identification of the black hole 
(BH) binary Cygnus X-1 (hereafter Cyg X-1) in the early 
1970's (e.g., Oda et al. 1971; Tananbaum et al. 1972; 
Thorne and Price 1975), X-ray observations have been 
playing an important role to reveal spectral and temporal 
properties of BH binaries, which are largely classified into 
two distinct states: the high/soft state and the low/hard 
state (e.g., Remillard and McClintock 2006; Done et 
al. 2007). In contrast to the high/soft state character- 
ized by the dominant disk emission (Mitsuda et al. 1984; 
Makishima et al, 1986) from the standard disk (Shakura 
& Sunyaev 1973), the spectrum in the low/hard state is 
expressed by a powerlaw with a photon index of ~ 1.5 
with an exponential cutoff at ~ 100 keV (e.g., Sunyaev & 
Triimper 1979) from a hot "corona" (e.g., Ichimaru 1977; 
Narayan & Yi 1995; chapter 8 in Kato, Fukue, and Mi- 
neshige 2008). Rapid time variabilities on a time scale 
of ~ ms (e.g., Miyamoto et al. 1991) only seen in the 
low/hard state have been studied in many ways (e.g., 
Nowak et al. 1999; Poutanen 2001; Pottschmidt et al. 
2003; Uttely et al. 2011; Torii et al. 2011), though the 
origin is still missing piece of puzzle, presumably due to 
observational difficulties of realizing both high sensitivity 
and large effective area. 

A distinctive approach is "shot analysis" (Negoro et 
al. 1994; Negoro 1995) adopted for Cyg X-1 obtained 
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with Ginga. This method is the time-domain stack- 
ing analysis to obtain universal properties behind non- 
periodic variability. It is in time-domain analysis that 
we can combine spectral information in a straightforward 
way. They found three main features: (1) the intensity 
changes time symmetrically, (2) both of the rise and de- 
cay curves are well represented by the superpositions of 
two exponential functions with time constants of ~ 0.1 s 
and r^ 1.0 s, and (3) the spectral variation is, by contrast, 
time asymmetric in the sense that it gradually softens to- 
ward the peak and instantly hardens across the peak (see 
Figure 2 in Negoro et al. 1994). These properties are fur- 
ther investigated with RXTE (Focke et al. 2005). The 
time constant of ~ 1 s far exceeds the local (dynamical or 
thermal) timescale of the innermost region, and should 
thus reflect accreting motion of gas element. Manmoto 
et al. (1996) proposed an interesting explanation that 
inward- forwarding accreting blob, causing an increase in 
X-ray flux, are reflected as sonic wave when it reaches the 
BH (Kato, Fukue, and Mineshige 2008), though further 
observational constraints have been awaited. 

The extension of this approach towards higher higher 
energies, ~ 200 keV, should be crucial because it may 
provide a hint to the physics causing the rapid spectral 
variation. Thus, we observed Cyg X-1 in the low/hard 
state with Suzaku (Mitsuda et al. 2007), by utilizing both 
the XIS (Koyama et al. 2007) located on the focus of the 
X-ray mirror (Serlemitsos et al. 2006) and the Hard X- 
ray Detector (HXD: Takahashi et al. 2007; Kokubun et 
al. 2007; Yamada et al. 2011) (Takahashi et al. 2007; 
Kokubun et al. 2007; Yamada et al. 2011). The distance, 
the mass, and the inclination of Cyg X-1 are 1.86^0 ii 
kpc (Reid et al. 2011; Xiang et al. 2011), 14.8 ± l.OM©, 
and 27.1 ±0.8° (Orosz et al. 2011), respectively. It has an 
09.7 lab supergiant, HD 226868 (Gies & Bolton 1986) 
with an orbital period of 5.599829 days (Brocksopp et 
al. 1999). Unless otherwise stated, errors refer to 90% 
confidence limits. 
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Fig. 1. — (a) A small portion for 3 s of a 0.5-10 keV light curve of Cyg X-1 taken with XISO operated in the P-sum mode. The time-bin 
size is 0.1 s. The arrows indicate the peaks judged as the shots, (b) The XISO shot profile created from all the light curve. 



2. OBSERVATION AND DATA REDUCTION 

Cyg X-1 data taken with Suzaku on 2008 April 18 (Ob- 
sID=403065010) are used in this letter, which is one of 
the 25 observations in its low/hard state (see Yamada 
2011 in details). The XISO was operated in the timing 
mode, or the Parallel-Sum (P-sum) mode, which is one of 
the clocking modes in the XIS. A timing resolution of the 
P-sum mode is ~ 7.8 ms. Data reduction of the timing 
mode are different from the standard one in the Grade 
selection criteria^: Grade (single event). Grade 1, and 
2 (double events) are used. The XIS background is not 
subtracted because it is less than 0.01 % of the signal 
events. 

The HXD data consisting of the PIN (10-60 keV) and 
GSO (50-300 keV) events are processed in the same man- 
ner as Torii et al. (2011). The events are selected by 
the criteria of elevation angle > 5°, cutoff rigidity > 6 
GV, and 500 s after and 180 s before the South Atlantic 
Anomaly. The non X-ray background (NXB) of the HXD 
modeled by Fukazawa et al. (2009) are subtracted from 
the HXD data. The NXB model reproduces the black- 
sky data with an accuracy of ~ 1% when the exposure is 
longer than ~ 40 ks. The cosmic X-ray background can 
be ignored in our analysis, since its contribution is less 
than ^ 0.1%. The simultaneous exposure for the XISO, 
PIN and GSO is 33.9 ks. 

3. DATA ANALYSIS AND RESULTS 

3.1. Definition and preparation 

We formulate the definition of the "shot" . Here t and 
C{t) is the event arriving time and the count rate at t. 
T is an interval of time over which C(t) is averaged, and 
C{t)j. denotes the average count rate over an interval of 
t — T < t < t -\-T. ta and ^5 refer to any time after and 
before t, satisfying the conditions of t < ta < t -\-T and 
t — T < th < t^ and / is the dimensionless parameter of 
order unity for the threshold. The peak time tp of each 
"shot" is a local maximum in C{t) defined as. 



{tp I C{U) < C{t) > C{ta) & C{t) > fC{t)^} (1) 

Equation (1) is nearly the same as that used in Negoro 
et al. (1994) and Focke et al. (2005), which works ro- 



^ see |http://www. a stro. isas.ac.jp/suzaku/analysis/xis/l 
psum_recipe/Fsum-recipe-2UlUU724.pdt 



bustly since no iteration is incorporated in this process. 
The only caveat is that when accidentally more than 
two adjacent bins of C{t) have the same value at the 
peak, the first of them is selected as the peak; e.g., when 
C{ti) < C{t2) = C{ts) > C{t4) is realized, ^2 is to be 
the peak. It is possible to impose optionally another 
constraint on the separation of time between the two 
successive peaks At; e.g.. At > gT works to avoid ac- 
cumulating the small peaks or fake events caused by the 
Poisson statistics, where g is the dimensionless param- 
eter of order unity. Since we aimed at accumulating as 
many photons as possible to quantify spectral change, 
we adopted ^ = 1, which means that a minimum of At 
equals T. 

An example of the light curve segment of XISO is shown 
in Figure 1(a). The events with 7.8 ms time resolution 
were binned into 0.1 s bins, resulting in 20-60 counts per 
bin. Intensity changes by a factor of ~ 2, so that Poisson 
fluctuation (~ 15% at la) is sufficiently smaller than the 
intrinsic variation. Employing / = 1.0 and T = 1.5 s, we 
actually applied the procedure of Equation (1) to the en- 
tire P-sum light curve, and identified 7524 shots in total. 
The distribution of At becomes a grossly exponential dis- 
tribution in agreement with the previous reports (Focke 
et al. 2005). Figure 1(b) shows the shot profile obtained 
by stacking the 7524 shots with reference to each peak. 
We can see two exponential slopes in the shot profile; 
shorter and longer decay time constants of ~ 0.1 s and 
1-2 s. Our primary focus is to extend this analysis to the 
HXD band, so we do not further investigate the shot pro- 
file and the XISO spectra due to incomplete calibration 
of the P-sum mode. 

3.2. Shot profiles of the HXD data 

According to the peak time determined with the XISO 
light curve, we have accumulated the PIN and GSO 
events and their NXB events^. To estimate pileup effects 
in a phenomenological way, we purposely tried to use ei- 
ther V inside or outside image core of XISO to obtain the 
shot profiles based on Yamada et al. (2012), though the 
shot profiles were not significantly changed. To investi- 
gate the energy dependence of the shots, we have utilized 
four energy bands: 10-20 keV and 20-60 keV from PIN, 

^ http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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and 50-100 keV and 100-200 keV from GSO. After sub- 
tracting the NXB events from the data and correcting 
them for dead time, we have obtained the 10-200 keV 
shot profiles with the HXD data. The stacked shot pro- 
files are divided by the count rates averaged over -4 to 
-2 s and 2 to 4 s to approximately correct them for the 
differences in the efficiency. 

The normalized shot profiles are shown in Figure 2(a)- 
(d). The derived profiles appear all very similar in shape 
to the one of XISO. However, energy dependencies are 
certainly found when their widths of the peaks are care- 
fully inspected. In Figure 2(a) and 2(d), the peak value 
in 10-20 keV is - 1.7 while - 1.5 in 100-200 keV, in- 
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Fig. 2. — The stacked profiles of the shots in the HXD bands by 
using the XISO light curve as a reference, (a)-(d) The background- 
subtracted and deadtime-corrected shot profiles in 10-20 keV, 20- 
60 keV, 50-100 keV, and 100-200 keV, which are renormalized by 
the individual averages over -4 to -2 and 2 to 4 s. (e)-(f) The 
hardness ratios of the shot profiles. The profiles in (b)-(d) are 
divided by the 10-20 keV profile in (a). 



dicating that the general trend that the higher photon 
energy is, the lower becomes the peak. This is neither 
due to incorrect background subtraction nor decrease in 
the sensitivity of the HXD, because the systematic un- 
certainty in the NXB subtraction is at most ~ 3% of 
the signal intensity even in the 100-200 keV band, and 
because we are referring to relative changes, instead of 
absolute values. To clarify the differences among these 
profiles, we divided the normalized shot profiles in the 
higher three bands by that in the 10-20 keV. As shown 
in Figure 2(e)-(g), the hardness ratios (relative to 10-20 
keV) gradually decrease towards the peak, but suddenly 
return to their average values immediately after (within 
0.1 s) the peak. Although this feature has been found in 
energies below ~ 60 keV in Negoro et al. (1994), we have 
not only confirmed the same trend up to ~ 200 keV, but 
also found that the spectral change is more prominent in 
higher energy of E > 100 keV. 

3.3. Quantification of the shot-phase-resolved spectra 

We then quantified its spectral change by accumulat- 
ing the HXD spectra according to the shot phase. The 
NXB events were accumulated in the same ways and sub- 
tracted. Figure 3 shows three examples of the derived 
shot-phase-resolved HXD spectra, corresponding to 0.15 
s before, right on, and 0.15 s after the peak. The expo- 
sure at the peak is 752.4 s (7524 shots x 0.1s). To grasp 
their characteristics in a model-independent way, we su- 
perposed the time-averaged spectrum, and show the ratio 
of the shot spectra to it in Figure 3. Aa shown evidently 
in Figure 3(b) by a clear turnover of the ratio above ~ 
100 keV, a spectral cutoff at the peak is lower than the 
averaged one. Furthermore, the spectral ratio before the 
peak shown in Figure 3(a) appears downward, while that 
after the peak is almost fiat, which is consistent with the 
gradual softening before the peak and instant hardening 
at the peak as seen in Figure 2. 

To consider physics underlying this spectral evolution 
we have fitted the 13 shot-phase-resolved HXD spectra 
with a typical model of Comptonization, compps (Pouta- 
nen & Svensson 1996) in the same manner as that in Torii 
et al. (2011). The seed photon is assumed to be a disk 
black body emission (Mitsuda et al. 1984; Makishima et 
al. 1986; Makishima et al. 2008) with a temperature of 
0.2 keV. The free parameters in the fits are the electron 
temperature Te, the optical depth r or the Compton y 
parameter, and the normalization A/'dbb- Note that if r 
is fixed, Tq is affected more by a spectral slope than a 
spectral cutoff. To avoid such a misunderstanding, we 
kept both r and Tq left free. As the shot-phase-resolved 
spectra do not have sufficient photon statistics, we fixed 
the reflection fraction Q at the value of 0.235, because 
the obtained value from the time-averaged spectrum is 
0.235^0.020- This implies that we assumed that the re- 
flection follows the primary continuum within ~ 0.1 s. 
The fits to all the spectra have been successful, result- 
ing in the best-fit parameters in Table 1. Even when 
considering the systematic error of the NXB in the GSO 
spectra, its contributions to the resultant values are less 
than - 1%. 

As the count rate increases on ~ s time scale, Tq and y 
decrease while r increases; when the count rate starts to 
decrease, all the parameters appear to return to the aver- 
aged values. To visualize this, we plot in Figure 4 the de- 
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Fig. 3. — The background-subtracted uFu spectra of the HXD, accumulated over different shot phases (black). The time- aver aged 
spectrum is given in red. Panel (a), (b), and (c) show the spectra integrated from -0.25 to -0.05 s before the peak, from -0.05 to 0.05 s 
around the peak, and from 0.05 to 0.25 s after the peak, respectively. Lower panels show the ratios to the time-averaged spectrum. 




Fig. 4. — Time evolution of the parameters of the Comptoniza- 
tion, determined by fitting the 12-300 keV HXD spectra in the 13 
shot phases with the Comptonization model. From top to bottom, 
the electron temperature, the optical depth, and the Compton y 
parameter are presented. The 90% range specified by the time- 
averaged spectrum is superposed by dotted and dashed lines. 



trary than the stacking analysis, but phase information 
is lost in the FT analysis. Further FT methods require 
more photons than a stacking method. Thus, we chose to 
used the stacking method and successfully extended the 
higher energy limit of the shot analysis up to ~ 200 keV, 
by utilizing the HXD data as well as the P-sum mode 
of the XIS. What we found are summarized as follows: 
(1) the shot feature is found at least up to ~ 200 keV 
with high statistical significance, (2) the shot profiles are 
approximately symmetric, though the hardness changes 
progressively more asymmetrically toward higher ener- 
gies of £; > 100 keV, and (3) the 10-200 keV spectrum 
at the peak shows lower energy cutoff than the time- 
averaged spectrum. By quantifing this feature in terms 
of the single-zone Comptonization, we found that (4) as 
a shot develops toward the peak, y and Tq decrease, while 
r and the fiux increases, and immediately past the shot 
peak, Tq and r (and hence y) suddenly return to their 



TABLE 1 

The fitting results of the shot-phase-resolved spectra. 



rived parameters in Table 1, as well as the time- averaged 
ones. Since our composite shot profile comprises a large 
number of relatively small individual shots, the averaged 
parameters are close to those at ~ 1 s from the peak. 
The gradual decrease in the ^/-parameter before the peak 
is consistent with the hardness decrease as seen in Fig- 
ure 2. The decrease in Tq around the peak clearly refiect 
the trend that the high-energy cutoff appears lowered at 
the peak as seen in Figure 3(b). Thus, the fitting re- 
sults are consistent with the hardness ratios in Figure 2 
and the spectral ratios in Figure 3. Note that A^dbb also 
increases along the shot profile, though we could not con- 
fidently measure the inner radius without using the soft 
X-ray data (cf. Makishima et al. 2008). 

4. DISCUSSION AND SUMMARY 

We performed the shot analysis to extract important 
information on understanding rapid hard X-ray variabil- 
ity, which can not be obtained by the Fourier transform 
(FT) methods (cf. Negoro et a. 2001, Legg et al. 2012, 
Torii et al. 2011). In general, FT methods are less arbi- 
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^ In a unit of 10^ R?^D^q cos 0, where R[n, ^10, and are the 
radius (km), the distance (10 kpc), and the inclination. 

^ The degree of freedom is 83 for the phased-sorted spectra, 
and 129 for the time-averaged one. 



RAPID SPECTRAL CHANGE IN HARD X-RAY OF CYGNUS X-1 WITH SUZAKU 



time-averaged values. 

Let us consider a possible physical mechanism to ex- 
plain the new findings and the previously- known features 
as well. The shot profile does not show any plateau at 
least down to ~ ms (Focke et al. 2005), which means 
that most of the luminosity is released almost time- 
symetrically within ~ 1 s or much shorter. Meanwhile, 
the hardness changes instantly within 0.1 s as shown 
in Figure 2, or much shorter as shown in Negoro et 
al. (1994). These features could be explained by some 
physical impulse or a discrete phenomenon, which can 
change properties of the radiation source in a short time. 

When accreting matter is assumed to be an ideal and 
non-relativistic gas, the entropy of the accreting gas, 5, 
with a temperature T and the density p is proportional 
to ln(P//)^) = ln(T//)^~-^), where 7 is the ratio of spe- 
cific heat capacities (5/3 for monatomic gas). It can be 
interpreted that the entropy decreases in some way as 
the fiux increases, but instantly increases at the peak, 
and returns to the mean value after the peak. This sug- 
gests the existence of some instant mechanism for direct 
entropy production (or heating). 

One of the possible ideas on the rapid intensity change 
have been considered as magnetic fiares analogous to the 
solar corona (Galeev et al. 1979), and recently more so- 
phisticated (cf., Poutanen & Fabian 1999; Zycki 2002). 
The magnetic fields are amplified by the differential ro- 
tation of the disk, and rise up into the corona where they 
reconnect and finally liberate their energy in fiare, mak- 
ing electrons accelerate. A typical magnetic model as- 
sumes so-called "avalanche magnetic fiare" , in which each 
fiare has a certain probability of triggering a neighboring 
one, producing long avalanches (Mineshige et al. 1994; 
Lyubarski, et al. 1997). An election temperature is ex- 
pected to increase if magnetic reconnection accelerate 
electrons or protons. It contradicts with the gradual 
decrease in Tq before the peak, while agrees with the 
instant increase at the peak and an expected time scale 
of ~ 1-100 ms when the energy dissipation occurs within 
r < 50i?g (i?g is the gravitational radius). However, it is 
still unknown how much magnetic power is stored in the 
corona. Magnetic reconnection in a plunge region may be 
likely based on the global three-dimensional MHD simu- 



lation (Machida and Matsumoto 2003). 

The mass propagating model (Manmoto et al. 1996, 
Negoro et al. 1995) gives an alternative explanation, be- 
cause the viscous time scale of the corona is ~ 1 s at 
100 i?g, and the only ~ 20% initial perturbation of mass 
accretion rate at ~ lOOi^g can change the luminosity 
by ~ 60%. Furthermore, the mass accretion refiected 
as a sonic ware can create the latter half of the peak; 
in the theoretical view point, the fiow passing a Bondi- 
type (not a disk-type) critical point does not always fall 
to the BH due to the strong outward centrifugal force 
unless angular momentum is very small (Kato, Fukue, 
and Mineshige 2008). The perturbation would start from 
overlapping region between the disk and corona, presum- 
ably ~ 50-100 Rg (Makishima et al. 2008), where intense 
turbulence is expected due to a large gap of the pressures 
and temperatures. Since the surface density of the corona 
is about four orders of the magnitude smaller than that of 
the disk, a little mass transfer from the disk to the corona 
in the overlapping region could increase mass accretion 
into the coronae. Shock phenomena might be possibly 
working for some reason because there are two or more 
sonic points around a BH (Nagakura and Yamada 2009). 

A rapid heating, such as magnetic reconnection, can 
explain the short (~ 0.1s) timescale of the shots, though 
the long (~ Is) timescale would be related to mass accre- 
tion time scale. Further observational studies are needed 
to completely understand the physics causing the rapid 
variability. For instance, shot profiles with distinct fea- 
tures, such as polarization (Laurent et al. 2011) or a 7- 
ray emission (Ling et al. 1987), could provide a new hint, 
which would be precisely measured by new missions like 
GEMS (Black et al. 2010) and ASTRO-H SGD (Taka- 
hashi et al. 2010). 
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